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1. Introduction 

The QCD radiation from incoming partons forms an inescapable component of the final 
state in all hard scattering processes at hadron colliders. This radiation leads to hadron 
formation that complicates the interpretation of events in a number of ways: by generating 
extra jets, by contaminating other jets, by modifying event shapes and global observables, 
and by changing the distributions of the products of the hard process. This last effect 
has been studied in great detail for the processes of electroweak boson production, with 
the result that the transverse momentum and rapidity distributions of W, Z and Higgs 
bosons at the Tevatron and LHC are predicted with good precision.^ The predictions for 

^See [1-3] and references therein. 



- 1 - 



the transverse momentum (qt) distributions in particular include resummation of terms 
enhanced at small qt to all orders in ag, matched with fixed-order calculations at higher 
qt values. The transverse momentum of the boson arises (neglecting the small intrinsic 
transverse momenta of the partons in the colliding hadrons) from its recoil against the 
transverse momenta of the radiated partons: qt = where 

Qt = -^Pti- (1-1) 

i 

The resummation of enhanced terms therefore requires a sum over emissions i subject to 
the constraint ( |1.1[ ) , which is most conveniently carried out in the transverse space of the 
impact parameter b Fourier conjugate to qrp: 

One then finds that the cumulative distribution in 6 = |b| contains terms of the form 
a^ln^{Qb), where Q is the scale of the hard process, set in this case by the mass of 
the electroweak boson, and p < 2n. These terms, which spoil the convergence of the 
perturbation series at large b, corresponding to small qt, are found to exponentiate [4-9]: 
that is, they can be assembled into an exponential function of terms that are limited to 
p < n + 1. This resummation procedure improves the convergence of the perturbation 
series at large values of b and hence allows one to extend predictions of the qt distribution 
to smaller values. 

Together with its vector transverse momentum prp^, every emission generates a con- 
tribution to the total hadronic transverse energy of the final state, Et, which, neglecting 
parton masses, is given by 

ET = Y,\PTi\- (1-3) 

i 

To first order in as (0 or 1 emissions) this quantity coincides with qx, but they differ in 
higher orders. In particular, at small qt there is the possibility of vectorial cancellation 
between the contributions of different emissions, whereas this cannot happen for the scalar 
Et- Thus the distribution of Et vanishes faster at the origin, and its peak is pushed 
to higher values. To resum these contributions at small Et, one should perform a one- 
dimensional Fourier transformation and work in terms of a 'transverse time' variable r 
conjugate to Et- 

5{ET-Y.\PTi\) = ^ [ dre-^-^n^^'"'"''- (l-^) 

Since the matrix elements involved are the same, one finds a similar pattern of enhanced 
terms at large r as was the case for large b: terms of the form ag In^(Qr) with p < 2n, which 
arise from an exponential function of terms with p < n + 1. Evaluation of the exponent to 
a certain level of precision (leading-logarithmic, LL, for p = n + 1, next-to-leading, NLL, 
for p = n, etc.) resums a corresponding class of enhanced terms and extends the validity 
of predictions to lower values of Et- 
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The resummation of Et in this way has received httle attention since the first papers 
on this topic, over 20 years ago [10-12]. This is surprising, as most of the effects of QCD 
radiation from incoming partons mentioned above depend on this variable rather than qx- 
A possible reason is that, unlike qt, Et also receives an important contribution from the 
so-called underlying event, which is thought to arise from secondary interactions between 
spectator partons. At present this can only be estimated from Monte Carlo simulations 
that include multiple parton interactions (MPI). Nevertheless it is worthwhile to predict as 
accurately as possible the component coming from the primary interaction, which carries 
important information about the hard process. For example, we expect the Et distributions 
in Higgs and vector boson production to be different, as they involve primarily gluon-gluon 
and quark-antiquark annihilation, respectively. Accurate estimates of the primary Et 
distribution are also important for improving the modelling of the underlying event. 

In the present paper we extend the resummation of Et in vector boson production 
to next-to-leading order (NLO) in the resummed exponent, parton distributions and co- 
efficient functions, and present for the first time the corresponding predictions for Higgs 
boson production. In Section ^ the resummation procedure is reviewed and extended to 
NLO; results on the resummed component are presented in Sect. ^. This component alone 
is not expected to describe the region of larger Et values, of the order of the boson mass; in 
Sect. 1^ we describe and apply a simple procedure for including the unresummed component 
at order ag ■ Section ^ presents Et distributions generated using the parton shower Monte 
Carlo programs HERWIG [13] and Herwig++ [14], which are compared with the analytical 
results and used to estimate of the effects of hadronization and the underlying event. Our 
conclusions are summarized in Sect. ^. Appendix A gives mathematical details of a com- 
parison between the resummation of the transverse energy Et and transverse momentum 
qt and Appendix B shows results for the LHC at lower centre-of-mass energy. 



2. Resummation method 



2.1 General procedure 



Here we generalize the results of ref. [11] to NLO resummation. The resummed component 
of the transverse energy distribution in the process /ii/i2 — )• FX at scale Q takes the form 



dap 



dQ^ dET 



— ^ / dxi dx2 dr e-'^^T fa/hAxi,^J-) fb/h2{x2,IJ') 

^ ah 



■ W[,ixiX2s;Q,T,f,) (2.1) 

where /a//i(a^,A*) is the parton distribution function (PDF) of parton a in hadron h at 
factorization scale ^, taken to be the same as the renormalization scale here. In what 
follows we use the MS renormalization scheme. As mentioned earlier, to take into account 
the constraint that the transverse energies of emitted partons should sum to Et, the 
resummation procedure is carried out in the domain that is Fourier conjugate to Et, using 
Eq. ( |1.4D . The transverse energy distribution (|2.1| ) is thus obtained by performing the 
inverse Fourier transformation with respect to the transverse time, r. The factor is 
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the perturbative and process-dependent partonic cross section that embodies the all-order 
resummation of the large logarithms ln((5r). Since r is conjugate to Et, the limit Et <C Q 
corresponds to Qt » 1. 

As in the case of transverse momentum resummation [15], the resummed partonic cross 
section can be written in the following universal form: 

Wab{s;Q,r,fJ,) = dzi / dz2Ccaias{^J'),Zl■,T,^^)Ccbias{^J■),Z2■,T,fJ,)6{Q'^-ZlZ2s) 

^ Jo Jo 

■ aUQ,asmSc{Q,r) . (2.2) 



Here cr^ is the cross section for the partonic subprocess c + c — )• F, where c,c = q,q (the 
quark qf and the antiquark qf/ can possibly have different flavours /, /') or c, c = g,g. The 
term Sc{Q, t) is the quark (c = q) or gluon (c = g) Sudakov form factor. In the case of Et 
resummation, this takes the form [11,12] 

5e(Q,r)=exp|-2^'^^ 2A,{as{q)) In^ + B,{asiq)) (l-e'"")}, (2.3) 

with c = q 01 g. The functions Ac{as), Bc{as), as well as the coefficient functions Cab in 
Eq. ( ^.2[ ), contain no In(Qr) terms and are perturbatively computable as power expansions 
with constant coefficients: 



oo 

'as 



n=l 

oo 

n=l 

oo 

Cabias,z) = 6ab6{l- z) + ^[^y C^\z) . (2.6) 



n=l 



Thus a calculation to NLO in as involves the coefficients a'^\ aI?\ , b'^c^ and C^J^ 
All these quantities are known for both the quark and gluon form factors and associated 
coefficient functions. Knowledge of the coefficients A^^^ leads to the resummation of the 
leading logarithmic (LL) contributions at small Et, which in the differential distribution 
are of the form cto^ln^ {Q / Et) / Et where p = 2n — 1. The coefficients i?*-^^ give the next- 
to-leading logarithmic (NLL) terms with p = 2n — 2, A^'^^ and C^^^ give the next-to-next- 
to-leading logarithmic (N^LL) terms with p = 2n — 3, and B^"^^ gives the N^LL terms with 
p = 2n — 4. With knowledge of all these terms, the first term neglected in the resummed 
part of the distribution is of order ag ln{Q / Et) / Et ■ 

In general the coefficient functions in Eq. (^]^) contain logarithms of jiT, which are 
eliminated by a suitable choice of factorization scale. To find the optimal factorization 
scale, we note that, to NLL accuracy, 

^lnPg(l-e'''0 ^ r ^InPq, (2.7) 

Q Jiro/r Q 
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where tq = exp(— 7^;) = 0.56146..., "fE being the Euler-Mascheroni constant. Therefore 
the effective lower hmit of the soft resummation is itq/t, and the parton distributions 
and coefficient functions should be evaluated at this scale. However, evaluation of parton 
distribution functions at an imaginary scale using the standard parametrizations is not 
feasible. We avoid this by noting that 



fa/h{x,q) =^ I —Uab{z\q,q)fh/h{x/z,q) {21 

I. J X Z 



b 

where Uab is the DGLAP evolution operator. Therefore 



/■I dz 

fa/hix,itJ.)= —Uabiz;ilJ;f^)fb/hix/z,ll) (2.9) 

Jx Z 

where the evolution operator Uab{z;ifJ., fi) is given to NLO by 

Uabiz; i/u, fi) = 6ab + -as(At) Pab{z) , (2.10) 

Pabiz) being the leading-order DGLAP splitting function. Similarly, in the coefficient 
functions we can write as{ifi) in terms of as(/^) using the definition of the running coupling: 

where (3 (as) = —ba^ + 0{a^), so that 

asiifi) = asifi) - iTrb[as{fi)]^ + 0(4) . (2.12) 



Furthermore, as the expressions ( p.l| ) and (|2.2| ) are convolutions, we can transfer the extra 
terms from ( |2.9D into the coefficient functions to obtain 



W^,{s;Q,t) = [\z, [\ 
^ Jo Jo 



dZ2 Cca{asiTo/T),Zi) Cg6(as (tq/t) , Z2) 5(Q^ - ZIZ2S) 

Jo Jo 

■ aUQ,asmSc{Q,r) (2.13) 

where 



Ccaiasif^),z) =Y —Ccdias{if^),z/z')Udaiz';itJ-,n). (2.14) 
d ^ 

Now the lowest-order coefficient function is of the form 

C(^{z) = C(PJ{z) = S,J{l-z) (2.15) 

and therefore 

CgHz) = CgHz) + i^P,.{z) . (2.16) 
Putting everything together, we have 

+ 00 



dap 



1 r+co 

— ^ j dr e-'^^^S,{Q,T) R,{s-Q,t) a!^,{Q,asm (2-17) 
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where, taking all PDFs and coefficient functions to be evaluated at scale /i = tq/t, 



i?,(s;g,r) = V / 

a,b 



dxi dx2 dzi 

Xl X2 Zl 



fa/h^{xi) fb/h2{x2) 



Cca{zi)C-,J \ . (2.18) 

\Z1X1X2S J 



To write as an integral over r > only, we note from (|2.9| ) and ( ^.lOl ) that when 

r — )• — r, to NLO the real parts of fa/hi and fb/h2 unchanged but the imaginary parts 
change sign. All other changes in ( ^.181 ) are beyond NLO. Thus, writing 



Rc^^ is symmetric with respect to r and R^^^ is antisymmetric. Defining 



(2.19) 



we therefore obtain 
dap 



21 ^ 

Q 

21 ^ 



dQ"^ dEr 



TTS 



2A,{as{q)) + B,{as{q)) 

Q 

2A,ias{q)) ln^ + Se(as(g)) 



V / dr e-^^*"" (s; Q, r) cosjF]^) (Q, r) - tEt} 

c -^0 ^ 

-Ri'\s;Q,T)sm{FP{Q,T)-TET}] af,(Q, as(Q)) (2.21) 
where, inserting ( |2.15| ) and ( p.l6| ) in ( |2.18| ) and defining ^ = Q'^/s, we have to NLO 



(1 — cos qt) 



sm gr 



(2.20) 



4«)(s;g,T)=4«)(e = QVs,T 
dxi dx2 



Xl X2 
+fc/h^{xi)falh2{x2)C, 



X1X2 



^XiX2 

Ri'\s;Q,r) = RPiC = QVs,r 
as sr-^ f dxi dx2 ~ 



Xl X2 
+fc/hi(.Xi)fa/h2{x2)Pc 



fa/hA^^)fc/h2{x2)Pa 



X1X2 



X1X2 



(2.22) 



It will be more useful to write, for example 

dxi dx2 



Xl X2 

dxi dx2 

Xl X2 

dxi dz 



fa/hi{xi)fc/h2{x2)Pca 



dz6 { z 



X1X2 



X1X2, 

fa/hi{xi)fc/h2{x2)Pca(.z) 



= [ — — fa/hA^l)fc/h2 (—] Pca{z) ■ 
J Xl Z ' ^ ' \ZXiJ 



(2.23) 
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This makes it more straightforward to interpret the +-prescription, which appears in some 
sphtting functions, as 

dxi dz 



I ^^^/'^^"'^i [^^'^"^ {-ill -^'^"^ a 

I' 'i^f'^/'^^i-^) [^/e-/.. (I;) {-I;) 



^ dxi [ i 



P{z) 
P{z) 



dzP{z) 



(2.24) 



We show in Appendix A that the results of resummation of the scalar transverse energy 
are identical to those of the more familiar resummation of vector transverse momentum at 
order as, as they should be since at most one parton is emitted at this order. 

The transverse energy computed here is the resummed component of hadronic initial- 
state radiation integrated over the full range of pseudorapidities rj. In ref. [11] the Et 
distribution of radiation emitted in a restricted rapidity range \r]\ < rj^ax was also es- 



timated. This was done by replacing the lower limit of integration in Eqs. ( 2.20 ) by 
Qc = <3 exp(— T/max), 1-6. assuming that radiation at q < Qc does not enter the detected 
region. This is justified at the leading-logarithmic level, where q/Q 6 ^ exp{—ri) and 
the scale dependence of the parton distributions and coefficient functions in Eq. ( p.lSf ) can 
be neglected. Then when ?7max = the form factor Sc is replaced by unity and Eq. (2.17) 
correctly predicts a delta-function at = times the Born cross section. However, this 
simple prescription cannot be correct at the NLO level, where the r dependence of the 
scale must be taken into account. Therefore we do not consider the Et distribution in a 
restricted rapidity range in the present paper. 

2.2 Vector boson production 

One of the best studied examples of resummation is in vector boson production through 
the partonic subprocess q + q' ^ V {V = W or Z): 



V 
'qq' 



(2.25) 



where at lowest order 



w 

qq' 



lV2GFM^\Vgg^ 



(2.26) 



with Vqqi the appropriate CKM matrix element and Vq,Aq the vector and axial couplings 
to the Z^. The coefficients in the quark form factor Sq{Q,T) are [8,16]: 



1 



vr 



9 2 
2 3__ 
16 



^(1) 



-Cf , 



(2.27) 



3C3 ) +CfCa\—7^~ 



11 



193 3^\ ^ fl7 TT^' 
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where is the Riemann C-function (C3 = 1.202 . . .), Cp = 4/3, Ca = 3, is the number 
of hght flavours, and 



6 / 9 ^ 



(2) 

The above expression for Bq is in a scheme where the subprocess cross section is given 
by the leading-order expression (|]2|). In the same scheme the NLO coefficient functions 
are [16, 17] 

Cgq{as,z) = {l + gCf(vr2-8)}5(l-z) + gC7^(l-z) 



(2.28) 



vr 



Cqg{as,z) 



as 
27r 



z{l-z) 



where the second line defines c'q \ The corresponding splitting functions are 



Cf 



— + -5{l - z) 

(1-2;)+ 2 



Pq,{z) = i + (1 _ . 

Equations (|2.22| )-( ^^ ) therefore give 

fq/hi{xi)fqlh. 



as 



1 (fa 
+ |/g/hl(2;i)/g7h2 

as /"""^ (ixi /"""^ dz 
XI Jo z 



Cf{1 - z) 



2 

+ 

OS 
2 



+ fq/hi{xi)fg/h 
\ZXl 



\ZXi 



2^(1-^) 



'^fq/hi{xi)fq/h 
fg/hM)fq/h2 + fq/hA^l)f9/h2 

|2CF/q/hi(a;i)/g/^ 

2CFfq/hAxi) ^fq/h^ Y 



^ fix 



+ S fg/hrixi)fg/h2 



ZXl 



+ fq/hi{Xl)fg/h2 



z 

ZXl 



2z 



(2.29) 



(2.30) 



(2.31) 
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2.3 Higgs boson production 

In the case of Higgs boson production the corresponding LO partonic subprocess is gluon 
fusion, g + g ^ H , through a massive-quark loop: 

aS(Q, asiQ)) = 5cg5-cg6{Q^ - ml)a^ , (2.32) 

where in the hmit of infinite quark mass 

28871^/2 

The coefficients in the gluon form factor Sg{Q,T) are [18-20] 

fif = (I + ii.^ - ^Ca) + ic. - 4) - fee. (2.34) 

(2) 

Here again, the above expression for Bg is in a scheme where the Higgs subprocess 
cross section is given by the leading-order expression (|]3|). In the same scheme the NLO 
coefficient functions are [21] 



Cgg{as,z) 



1 + 



as 
4tt 



vr 



Ca ( 2 - — 1 + 5 + 47r2 



5{1 - z) 



l + ^c«)<5(l-.) 



as . 



Cgq{a^,z) = Cgq{as,z) = —Cpz . 



27r 



The corresponding splitting functions are 

z 1 — z 



Pgg{z) = 2Ca 



+ 



(1 - z)+ 

l + {l-zf 



+ z{l-z) 



+ -{llCA-2nf)5{l- z) 





Pgq{z) = Pgq{z) = Cf 

z 

Equations ( |2.22| )-( ^3^ ) therefore give 



Xl 



1 + 



^2c(i) 



vr 



as 



TT 



dz 



+— — 



fg/hAxi)fs/h2 + fs/hM)fg/h2 



Icfz 



as dxi dz 



i Jo 



P99i^) 



+ 



fg/hAxi)fs/h, + fs/hr{xi)fg/h2 



(2.35) 



(2.36) 
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Figure 1: Resummed component of the transverse energy distribution in Z*^ boson production at 
the Tevatron and LHC. The curves show the effects of the coefficients in the quark form factor: 



black, all coefficients; magenta omitting Bq , blue Aq''' and A^f' only; red Aq''' only. 



l(2) 



2CAln( 1-^) +^inCA-2nf) 



dz 



z 1 — z 

H +z(l - z] 



1- z 



Cf 



l + {l-zf 



fg/hrixi)fs/h2 + fs/hiixi)fg/h2 

where fs = T^qUq + fg)- 

3. Resummed distributions 
3.1 Vector boson production 

Figure |l| shows the resummed component of the transverse energy distribution in boson 
production at the Tevatron (pp at ^/s = 1.96 TeV) and LHC {pp at ^/s = 14 TeV).^ For all 
calculations, we use the MSTW 2008 NLO parton distributions [22]. The different curves 
show the effects of the subleading coefficients ( p. 27 ) in the quark form factor. We see that 
while Bq^^ has a large effect (the difference between the blue and magenta curves), the 
effects of the other subleading coefficients are quite small. 

The peak of the resummed distribution lies at around Et ~ 10 GeV at the Tevatron, 
rising to ~ 14 GeV at the LHC. This is comfortably below Mz, justifying the resumma- 
tion of logarithms of Et/Mz in the peak region. However, at LHC energy the predicted 
distribution has a substantial tail at larger values of Et, indicating that the higher-order 
terms generated by the resummation formula remain significant even when the logarithms 
are not large. In addition, the LHC prediction does not go to zero as it should at small 
Et- However, this region is sensitive to the treatment of non-perturbative effects such as 



(2.37) 



^Results for pp at — 7 TeV are given in Appendix B. 
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Figure 2: Resummed component of the transverse energy distribution in Higgs boson production 
at the Tevatron and LHC. The curves shows the effects of the coefficients in the gluon form factor: 
black, all coefficients; magenta omitting Sg^"*; blue Ag^'^ and aI^^ only; red Ag^'^ only. 



the behaviour of the strong coupling at low scales (we freeze its value below 1 GeV) and 
the upper limit in the integral over transverse time (we set Tmax = 1/A where A is the 
two-loop QCD scale parameter, set to 200 MeV here). 

The resummed component for boson production looks very similar, apart of course 
from the overall normalization, and therefore we do not show it here. Predictions with 
matching to fixed order will be presented in Section |^. 



3.2 Higgs boson production 

Figure ^ shows the resummed component of the transverse energy distribution in Higgs 
boson production at the Tevatron and LHC, for a Higgs mass of 115 GeV. The effects of 



subleading terms in the gluon form factor (2.34) are more marked than those of the quark 



form factor discussed above. The distribution peaks at large values of Et, around 40 GeV 
at the Tevatron, rising to ~ 50 GeV at the LHC. This is due to the larger colour charge 
of the gluon. However, together with the large effects of subleading terms, it does make 
the reliability of the resummed predictions more questionable. Also in contrast to the 
vector boson case, the suppression at low and high Et is if anything too great, resulting 
in negative values below 16 GeV and above 120 GeV at Tevatron energy. 



4. Matching to fixed order 

The resummed distributions presented above include only terms that are logarithmically 
enhanced at small Ej-- To extend the predictions to larger Et we must match the resum- 
mation to fixed-order calculations. To avoid double counting of the resummed terms, the 
corresponding contribution must be subtracted from the fixed-order result. 

We consider here only matching to first order in as- To this order the Et distribution 
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Figure 3: Left: order-as Et distribution in Z" production at the Tevatron; solid, full prediction; 
dashed, fit to enhanced terms. Right: difference between full prediction and fit to enhanced terms. 




Figure 4: Predicted Et distribution in production at the Tevatron and LHC. Solid: resummed 
prediction matched to 0{as)- Dashed: resummed only. Dot-dashed: ©(as) only. 



for Et > has the form 



'^'^ -^iAlnET + B) + C{Et) (4.1) 



(IEt Et 



where A and B are constants (for a given process and collision energy) and the function 
C{Et) is regular at Et = 0. The terms involving A and B are already included in the 
resummed prediction, and therefore we have only to add the regular function C to it to 
obtain a prediction that is matched to the 0{a^) result. This function is determined by 
fitting the 0{as) prediction for ETda/dET to a linear function of \b.Et at small Et, 
extracting the coefficients A and i?, and then subtracting the enhanced terms in Eq. ([4.1|). 

4.1 Vector boson production 

The above matching procedure is illustrated for production at the Tevatron in Fig. |3|. 
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20 40 60 20 40 



Et (GeV) Et (GeV) 

Figure 5: Predicted Et distribution in W++W^ production at the Tevatron and LHC. Solid: 
resummed prediction matched to ©(as)- Dashed: resummed only. Dot-dashed: 0{as) only. 



The fit to the logarithmically enhanced terms gives excellent agreement with the order- 
as result out to around 20 GeV, confirming the dominance of such terms throughout the 
region of the peak in Fig. |l|. The remainder function C{Et) vanishes at small Et and rises 
to around 10 pb/GeV, falling off slowly at large Et- Consequently the matching correction 
to the resummed prediction is small and roughly constant throughout the region 40-100 
GeV, as shown in Fig. ^. 

As shown on the right in Fig. ^, the situation is similar at LHC energy: the matching 
correction is small, although in this case it is negative below about 40 GeV. The large tail 
at high Et and the bad behaviour at low Et, due to uncompensated higher-order terms 
generated by resummation, are not much affected by matching to this order. 

The corresponding matched predictions for boson production are shown in Fig. |5|. 
As remarked earlier, the form of the resummed distribution is very similar to that for 
boson production, and again the matching correction is small. 

4.2 Higgs boson production 

Adopting the same matching procedure for Higgs boson production, we find the results 
shown in Figs. |^ and |^. The form of the matching correction is similar to that for vector 
bosons, but its effect is rather different. The roughly constant, then slowly decreasing, 
correction in the region 20-100 GeV is not small compared to the resummed result and 
therefore it raises the whole distribution by a significant amount throughout this region. 
This has the beneficial effect of compensating the negative values at low and high Et at 
Tevatron energy. However, it further enhances the high-£'r tail of the distribution at LHC 
energy. This, together with the relatively large correction in the peak region, casts further 
doubt on the reliability of the predictions in the case of Higgs production. 

5. Monte Carlo comparisons 

In this section we compare the resummed and matched distributions obtained above with 
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Figure 6: Left: order-as Et distribution in Higgs boson production at the Tevatron; solid, full 
prediction; dashed, fit to enhanced terms. Right: difference between ftill prediction and fit to 
enhanced terms. 
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Figure 7: Predicted Et distribution in Higgs boson production at the Tevatron and LHC. Solid: 
resummed prediction matched to 0{as)- Dashed: resummed only. Dot-dashed: 0{as) only. 



the predictions of the parton shower Monte Carlo programs HERWIG [13] and Herwig++ [14]. 

Comparisons are performed first at the parton level, that is, after QCD showering 
from the incoming and outgoing partons of the hard subprocess. We say "incoming and 
outgoing" because both programs apply hard matrix element corrections: in addition to 
the Born process, order-as real emission hard subprocesses are included in phase-space 
regions not covered by showering from the Born process. 

After showering, the Monte Carlo programs apply a hadronization model to convert 
the partonic final state to a hadronic one. We show the effects of hadronization in the case 
of HERWIG only; those in Herwig++ are broadly similar since both programs use basically the 
same cluster hadronization model. The programs also model the underlying event, which 
arises from the interactions of spectator partons and makes a significant contribution to the 
hadronic transverse energy. In this case wc show only the underlying event prediction of 
Herwig++, since the default model used in HERWIG has been found to give an unsatisfactory 
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Figure 8: Predicted Et distribution in Z'^ boson production at the Tevatron and LHC. Comparison 
of resummed and Monte Carlo results. 

description of Tevatron data. For an improved simulation of the underlying event, HERWIG 
can be interfaced to the multiple interaction package JIMMY [23], which is similar to the 
model built into Herwig++. 

5.1 Vector boson production 

Figure ^ shows the comparisons for production at the Tevatron and LHC. The HERWIG 
predictions are renormalized by a factor of 1.3 to account for the increase in the cross 
section from LO to NLO. The Herwig++ results were not renormalized, because they were 
obtained using LO** parton distributions [24], which aim to reproduce the NLO cross 
section. We see that the parton-level Monte Carlo predictions of both programs agree fairly 
well with the matched resummed results above about 15 GeV, but Herwig++ generates a 
substantially higher number of events with low values of Et- A similar pattern is evident in 
the results on boson production, shown in Fig. The effects of hadronization, shown 
by the difference between the blue and magenta histograms, are also similar for both vector 
bosons. They generate a significant shift in the distribution, of around 10 GeV at Tevatron 
energy and 20 GeV at LHC. 

5.2 Higgs boson production 



As may be seen from Fig. 10, the agreement between the resummed and parton-level Monte 
Carlo results is less good in the case of Higgs boson production than it was for vector bosons. 
Here we have renormalized the HERWIG predictions by a factor of 2 to allow for the larger 
NLO correction to the cross section. Then the Monte Carlo Et distributions agree quite 
well with each other but fall well below the matched resummed predictions. Fair agreement 
above about 40 GeV can be achieved by adjusting the normalization, but then the Monte 
Carlos predict more events at lower Et- The effect of hadronization is similar to that in 
vector boson production, viz. a shift of about 10 GeV at the Tevatron rising to 20 GeV at 
the LHC, which actually brings the HERWIG distribution into somewhat better agreement 
with the resummed result. 
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Figure 9: Predicted Et distribution in W++W boson production at the Tevatron and LHC. 
Comparison of resummed and Monte Carlo results. 
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Figure 10: Predicted Et distribution in Higgs boson production at the Tevatron and LHC. Com- 
parison of resummed and Monte Carlo results. 



5.3 Modelling the underlying event 

Figures ^ and 12 show the parton-level Herwig++ predictions for the Et distribution 
in and Higgs boson production, respectively, with the contributions from initial-state 
radiation (in red, already shown in Figs. |^ and IC ), the underlying event (blue) and the 
combination of the two (green). The underlying event is modelled using multiple parton 
interactions; see ref. [14] for details. Clearly it has a very significant effect on the Et 
distribution. However, this effect is substantially independent of the hard subprocess, as 



may be seen from the comparison of different subprocesses in Fig. 13 



We find that the probability distribution of the Et contribution of the underlying 
event in the Herwig++ Monte Carlo can be represented quite well by a Fermi distribution: 




Figure 11: Predicted Et distribution in Z° boson production at the Tevatron and LHC. Monte 
Carlo results including underlying event. 
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Figure 12: Predicted Et distribution in Higgs boson production at the Tevatron and LHC. Monte 
Carlo results including underlying event. 



where the normalization is 



AT = Tin 



exp 



+ 1 



(5.2) 



The dependence of the "chemical potential" /i and "temperature" T on the hadronic colli- 
sion energy is shown in Fig. 14, The red curves show fits to the energy dependence of the 
form 



1 



where the coefficients in the fits are A = 20(1), B = 0.030(4), q = 36(2), r = 0.28(3). 



(5.3) 



6. Conclusions 



We have extended the resummation of the hadronic transverse energy Et in vector boson 
production to next-to- leading order (NLO) in the resummed exponent, parton distributions 
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Figure 13: Comparison of Et distributions of the underlying event in different subprocesses at 
the LHC. 
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Figure 14: Fitted values of the parameters of the underlying event in Higgs production in pp 
collisions at various energies. 

and coefficient functions, and also presented for the first time the corresponding predictions 
for Higgs boson production. We have matched the resummed results to the corresponding 
0{as) predictions, by adding the contributions in that order which are not included in the 
resummation. In addition we have compared with parton shower Monte Carlo results and 
illustrated the effects of hadronization and the underlying event. 

In the case of vector boson production, the resummation procedure appears stable 
and the parton-level results should be quite reliable. The leading-order mechanism of 
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quark-antiquark annihilation typically generates a moderate amount of transverse energy 
in initial-state QCD radiation. Consequently the effects of subleading resummed terms and 
fixed-order matching are small and the peak of the Ej- distribution lies well below the boson 
mass scale, where resummation makes good sense. The comparisons with Monte Carlo 
programs reveal some discrepancies but these are at the level of disagreements between 
different programs; in this case the resummed predictions should be more reliable (at 
parton level) than existing Monte Carlos. The programs suggest that the non-perturbative 
effects of hadronization and the underlying event are substantial. These effects can however 
be modelled in a process-independent way. We have suggested a simple parametrization of 
the contribution of the underlying event. 

The situation in Higgs boson production is not so good. The dominant mechanism of 
gluon fusion generates copious ISR and the effects of subleading terms and matching are 
large. The resummed Et distribution peaks at a value that is not parametrically smaller 
than the Higgs mass and the behaviour at low and high Et is unphysical before matching. 
The discrepancies between the matched resummed and Monte Carlo predictions are sub- 
stantially greater than those between different programs, even allowing for uncertainties in 
the overall cross section. All this suggests that there are significant higher-order corrections 
that are not taken into account, either further subleading logarithms or unenhanced terms 
beyond NLO. It would be interesting (but very challenging) to attempt to extract such 
terms from the available NNLO calculations of Higgs production. 
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A. Relation to transverse momentum resummation 



Here we demonstrate the equivalence of transverse energy and transverse momentum re- 
summation at order as- Expanding Eq. ( |2.3D to this order, using (2/7) and substituting 
into and ( |2.2D , we find terms involving the integrals 

MQ, Et) = ^ dr e"*^^- In^ (^) (A.l) 

with p = 1,2. At this order, evaluating the PDFs at the scale itq/t leads to single- 



logarithmic terms of the same form when we use (2^) to write 

^1 dz 

b 



fa/hix, iro/r)) = fa/dx, Q) - ^ In (^) ^ ^p^,(^)/,/,(^/^, Q) . (A.2) 
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The integral ( A.l ) may be evaluated from 



IpiQ,ET) = —I{Q,Et;u)\u=o 



where 



IiQ,ET;u) 



2lT 



Writing r = iz/ET, we have 

I{Q,Et;u) 



2ttEx KExTq 



dr e 



Q 



QZ_Y 



dz e^ . 



(A.3) 
(A.4) 

(A.5) 



We can safely deform the integration contour around the branch cut along the negative 
real axis to obtain 



Z{Q,Er,u) 



I ( Q Y 



TtEx \ExTq 

which, recalling that Inrg = —'Je = gives 



sin(7rn) r(l + u) , 



Q 



I,{Q,Et) = -^, UQ,Et) = -^ln[^ 



(A.6) 



(A.7) 



The resummed component of the transverse momentum {qx) distribution takes the 



form 



dap 



dQ^ dqr 



rl fl /"OO 

„ f, JO Jo Jo 

(Ai 



'■es- a,b 
■F 



■ WabixiX2s;Q,b) 
where 60 = 2exp(— 7^;), 



Wai,{s;Q,b) = S2 dzi / dz2 Cca{as{bo/b), zi) Ccb{as{bo/b), Z2) d{Q'^ - ziZ2s) 
Jo Jo 



and 



af,(Q,as(Q)) S,{Q,b) 



5e(Q,6) = exp<;-2 /"^ ^ 

ibo/b Q 



2A,iasiq)) In^ + B^Mq)) 



(A.9) 



(A.IO) 



Expanding to order as, we find the same terms as in the E^ resummation except that 
( |A.1|) is replaced by 



MQ^Qt) = QT / dbbJoibqx) InPiQb/bo) . 
Jo 

It therefore suffices to show that 

Xp{Q,qT) =Zp{Q,ET = qr) for p = 1, 2. 



(A.ll) 



(A.12) 
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Now corresponding to 



we have 



Using the result 



gives 



and hence 



I{Q,qT;u) = qr I dbbJoibgr) ( — 



' dtt^-' Mt) = ^ sin 



I{Q,qT;u) 



^i{Q,qT) 



/2Q_y 



1 



sm 



TTU 



\ p2 



2 / 



1 + 



, l2{Q,qT) = — —Ini — 

qr qx Kqt 



(A.13) 

(A.14) 
(A.15) 
(A.16) 



in agreement with ([A.?]) and ( A. 12 ). Notice, however, that the higher {p > 2) derivatives 
of I and I differ, corresponding to the difference between Ej- and qT resummation beyond 
0{as). 
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B. Results for the LHC at 7 TeV 




Figure 15: Predicted Et distributions in Z'^, W++W and Higgs boson production in pp collisions 
at = 7 TeV . 



We show here results for the LHC operating at a centre-of-mass energy of 7 TeV, corre- 
sponding to those shown earlier for 14 TeV. Apart from the normalization, the predictions 
for the two energies are very similar, with only a slight downward shift in the position of 
the peak in the Et distribution at the lower energy. 



- 22 - 



References 

[1] G. Bozzi, S. Catani, D. de Florian and M. Grazzini, Nucl. Phys. B 791 (2008) 1 

[arXiv:0705.3887 [hep-ph]]. 

[2] G. Bozzi, S. Catani, G. Ferrera, D. de Florian and M. Grazzini, Nucl. Phys. B 815, 174 
(2009) [arXiv:0812.2862 [hep-ph]]. 

[3] S. Mantry and F. Petriello, arXiv:0911.4135. 

[4] Y. L. Dokshitzer, D. Diakonov and S. I. Troian, Phys. Rep. 58 (1980) 269. 

[5] G. Parisi and R. Petronzio, Nucl. Phys. B154 (1979) 427. 

[6] G. Curci, M. Greco and Y. Srivastava, Nucl. Phys. B159 (1979) 451. 

[7] A. Bassetto, M. Ciafaloni and G. Marchesini, Nucl. Phys. B163 (1980) 477. 

[8] J. Kodaira and L. Trentadue, Phys. Lett. B112 (1982) 66, Phys. Lett. B123 (1983) 335. 

[9] J. C. Collins, D. E. Soper and G. Sterman, Nucl. Phys. B 250 (1985) 199. 

[10] F. Halzen, A. D. Martin, D. M. Scott and M. P. Tuite, Z. Phys. C 14 (1982) 351. 

[11] C. T. H. Davies and B. R. Webber, Z. Phys. C 24 (1984) 133. 

[12] G. AltarelU, G. Martinelli and F. Rapuano, Z. Phys. C 32 (1986) 369. 

[13] G. Corcella, L G. Knowles, G. Marchesini, S. Moretti, K. Odagiri, P. Richardson, 
M. H. Seymour and B. R. Webber, JHEP 0101 (2001) 010 [arXiv:hep-ph/0011363]; 
arXiv:hep-ph/0210213. http : //pro j ects . hepf orge . org/ f herwig/ 

[14] M. Bahr et al, Eur. Phys. J. C 58 (2008) 639 [arXiv:0803.0883 [hep-ph]]; arXiv:0812.0529 
[hep-ph] . http : // pro j ect s . hepf orge . org/herwig/ 

[15] S. Catani, D. de Florian and M. Grazzini, Nucl. Phys. B 596, 299 (2001) 
[arXiv:hep-ph/0008184]. 

[16] C. T. Davies and W. J. Stirling, Nucl. Phys. B244 (1984) 337. 

[17] C. Balazs, J. W. Qiu and C. P. Yuan, Phys. Lett. B 355 (1995) 548 [arXiv:hep-ph/9505203]. 

[18] D. do Florian and M. Grazzini, Phys. Rev. Lett. 85 (2000) 4678 [arXiv:hcp-ph/0008152]. 

[19] D. de Florian and M. Grazzini, Nucl. Phys. B 616 (2001) 247 [arXiv:hep-ph/0108273]. 

[20] S. Catani, E. D'Emilio and L. Trentadue, Phys. Lett. B211 (1988) 335. 

[21] R. P. Kauffnian, Phys. Rev. D45 (1992) 1512. 

[22] A. D. Martin, W. J. Stirling, R. S. Thorne and G. Watt, Eur. Phys. J. C 63 (2009) 189 
[arXiv:0901.0002 [hep-ph]]. 

[23] J. M. Butterworth, J. R. Forshaw and M. H. Seymour, Z. Phys. C 72 (1996) 637 
[arXiv:hep-ph/9601371]. http://projects.hepforge.org/jiinmy/ 

[24] A. Sherstnev and R. S. Thorne, Eur. Phys. J. C 55 (2008) 553 [arXiv:071 1.2473 [hep-ph]]; 
arXiv:0807.2132 [hep-ph]. 



-23- 



